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ABSTRACT 



We report angular diameters of HIP 560 and 21547, two F spectral type pre- 
main sequence members of the (3 Pic Moving Group. We used the East- West 
314-m long baseline of the CHARA Array. The measured limb-darkened angular 
diameters of HIP 560 and 21547 are 0.492 ± 0.032 and 0.518 ± 0.009 mas, 
respectively. The corresponding stellar radii are 2.1 and 1.6 R for HIP 560 and 
HIP 21547 respectively. These values indicate that the stars are truly young. 
Analyses using the evolutionary tracks calculated by Siess, Dufour, and Forestini 
and the tracks of the Yonsei-Yale group yield consistent results. Analyzing the 
measurements on an angular diameter vs color diagram we find that the ages of 
the two stars are indistinguishable; their average value is 13 ± 2 MY. The masses 
of HIP 560 and 21547 are 1.65 ±0.02 and 1.75 ±0.05 M , respectively. However, 
analysis of the stellar parameters on a Hertzsprung-Russell Diagram yields ages 
at least 5 MY older. Both stars are rapid rotators. The discrepancy between the 
two types of analyses has a natural explanation in gravitational darkening. Stellar 
oblateness, however, does not affect our measurements of angular diameters. 



Subject headings: stars: Pre-Main Sequence — stars: Beta Pic Moving Group 
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1. Introduction 

The age and mass of a pre-main sequence (PMS) star are usually estimated from its 
location in the Hertzsprung-Russell diagram (HRD) relative to theoretical calculations 
of stellar evolution. Unfortunately, differences among the theoretical calculations can 
produce mass and age estimates discrepant by factors of 2 to 3 (Hillenbrand and White, 
2004; Simon, 2006). The consequences are that the mass-spectrum of stars produced 
in a star-forming region, the region's star-forming history, and the chronology of planet 
formation are imprecisely known. The latter point has taken on a pressing importance as 
astronomers are poised to detect exoplanets associated with young stars. 

Measurement of the diameters of young stars as they contract to the main sequence 
can provide their ages and a new independent test of the theoretical PMS evolutionary 
tracks (Simon, 2006). Once a reliable set of evolutionary tracks is established, measurement 
of the diameters of PMS stars will yield measurements of their absolute ages. Remarkably, 
members of the (5 Pic Moving Group (BPMG) are young enough (10-20 MY) and near 
enough (< 50 pc; e.g. Torres et al., 2006) that some are resolvable with the optical/IR 
interferometric array of the Center for High Angular Resolution Astronomy (CHARA). 
The brighter members of the BPMG have parallaxes measured by HIPPARCOS. The ages 
of individual stars in the BPMG are determined by HRD fitting and the strength of the 
Li-line. An age for the group as a whole may be estimated dynamically by tracing the 
group's motion back to a common origin. The age estimates range from 11.2 ± 0.3 to 
22 ± 12 MY (see Table 2 in Fernadez et al. 2008; Yee and Jensen, 2010). Whether the 
spread represents uncertainties of the dating methods or a range in individual formation 
times is unknown. Astronomers already have evidence of planet formation associated with 
the members of the BPMG. Lagrange et al (2010) imaged the exoplanet of j3 Pic, an A-star. 
Debris disks are understood as indicators of planet formation (Quillen et al. 2007); f3 Pic 
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and AU Mic, a BPMG member of Ml spectral type, are known to have debris disks (Smith 
and Terrile 1984; Lagrange et al 2010; Liu 2004). 

We measured the angular diameters of two stars in the BPMG , HIP 560 and 21547, 
using the CHARA Array. §2 describes our observations and §3 presents our analysis and 
discussion of the diameter measurements. We summarize our results in 54. 



2. Interferometric Measurements of HIP 560 and HIP 21547 

We observed at the CHARA Array located on Mt. Wilson, CA, on the nights of UT 
11, 12, and 13 Sept. 2010. Our observing assignment was made through an allotment 
made available competitively to the astronomical community by CHARA and administered 
by the National Optical Astronomical Observatory through its observing time application 
process. The array consists of six 1-m telescopes in a Y-configuration on baselines from 34 
to 331 m (ten Brummelaar et al. 2005). We used the CLASSIC beam combiner operating 
in the H and K' bands to observe with the two telescopes on the 314 m El-Wl baseline. 
At the CHARA Array these bands are centered at 1.673 and 2.133 [im wavelengths, with 
bandwidths, 0.285 and 0.349 fim respectively. 

Table 1 lists the properties of HIP 560 and HIP 21547 we need for this paper; the 
parallactic distances are van Leeuwen et al.'s (1997) revised HIPPARCOS values. Neither 
star is known to be a spectroscopic binary. HIP 21547 is one component of a common 
proper motion binary; its companion, GJ 3305, is at about 1' angular separation. We 
checked the near-IR and thermal-IR fluxes of the stars as given by 2MASS and WISm 
Neither shows evidence of excess emision in the near-IR. HIP 560 is not included in the 
WISE Preliminary Data Release. For HIP 21547, the WISE magnitudes at 3.6, 12, and 22 
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fxm are consistent with the 2MASS K-band magnitude, ~ 4.5 (Table 1). HIP 21547 is ~ 
0.5 mag brighter in WISE band 2 at 4.5 fim; on a color-color diagram this still places the 
star in the region of normal stars without a debris disk. Both stars are rapid rotators with 
values of vsin i in the range typical of F stars (Abt and Hunter 1962). 

The essential observational datum of an interferometer measurement is the fringe 
contrast or visibility, V, of the target or calibrator. An observation of a program star or its 
calibrators took about 10 minutes each and consisted of a series of data scans through the 
central interferometer fringe and a sequence of shutters. Each program star observation was 
bracketed by observations of one or two stellar calibrators with angular diameters smaller 
than 0.3 mas and located within 10° of the target (Table 1). An observation of HIP 560 or 
21547 and their calibrators required about 2 hours to obtain good S/N. 

The calibrators were chosen to be unresolved by the interferometer, but their 
diameters are nonetheless finite (Table 1). S. Ridgway (priv. comm.) kindly estimated 
the limb-darkened diameters, <Pld, of the calibrators, using Kervella et al's (2004) angular 
diameter fits for main sequence stars as a function of their V-K color, K-band magnitude, 
the parallactic distance of the star, and an estimate of its extinction (always small, 
Ay < 0.07). We scaled the observed calibrator visibilities to those of a point-like star. We 
then calibrated the target visibilities by ratioing them to the scaled calibrator visibilities. 
Table 2 lists the wavelength, MJD time of observation, projected baseline in meters, its 
position angle on the sky, measured eastward from north, the calibrated visibilities V ca i 
of the targets, and their uncertainties avcai- The &vcai include the uncertainties of the 
measured visibilities of target and calibrators and an assumed ±10% uncertainty in the 
calculated diameters of the calibrators. 

Figures 1 and 2 show the calibrated visibilities at H and K vs spatial frequency 
(projected baseline/wavelength). Fig. 1 for HIP 560 indicates a systematic difference in 
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angular diameters measured at H and K. We do not think this can be attributed to a 
property of the star because HIP 21547, also an F spectral type star, does not show the 
effect. The difference is probably measurement error. At angular diameters as small as 
those of HIP 560 and HIP 21547 the curves for a uniform disk and limb-darkened disk 
are nearly indistinguishable^. Nonetheless, we did analyze the visibilities with respect to 
limb-darkened models following Hanbury-Brown et al.'s (1974) analysis. Their expression 
for the visibility of a limb-darkened star can be written 



V LD {x) = [— g— + y] [gC 1 - u >) v ud + cosx )\ 

where Vjjd is the visibility of a uniform stellar disk, Vjjd = 2 : ^^, in which J\ is the Bessel 
function of order 1, and x = with (j) the star's angular diameter, B the projected 
baseline during the scan, and A the wavelength of observation. We used values of the limb 
darkening parameter, u\, derived by Claret et al (1995) appropriate to stars that have 
effective temperatures of 1.5 to 1.7 M stars 10 to 20 MY old, u x = 0.24 and 0.20, at H 
and K respectively. We included a ±10% uncertainty in the values of u\ in the calculation 
of the uncertainties of the measured stellar diameters. 

We fit each calibrated visibility to a limb-darkened diameter, averaged the individual 
values, and present the results, <pDiam{ ma s), for HIP 560 and 21547 in Table 3. The 
uncertainties are standard deviations of the mean. The uncertainties are dominated by the 
scatter of the individual visibility measurements. Table 3 also lists ^Diam the "absolute" 
angular diameter, the value of <fi scaled to a common distance of 10 pc and the corresponding 
stellar radii. The uncertainties in §Diam include the uncertainties in the HIPPARCOS 



2 This would not be true for stars more strongly limb-darkened than F-stars (van Belle et 
al. 2001). 



-7- 



distances (van Leeuven et al. 2007). 



3. Analysis and Discussion 
3.1. Analysis in the $ vs (V — K) and the HR Diagrams 



Figure 3 shows $r>?.am compared with diameters predicted by PMS evolution models 
calculated by Siess et al (2000) (SDF) and the Yonsei-Yale (Y2) models calculated by Yi et 
al (2003). The Y2 calculations^ provide the luminosity, L, and effective temperature, T e // 
at 5 MY intervals during the contraction; we calculated the corresponding radii through 
the defining relation L = 4nR 2 aT^jr where a is the Stefan-Boltzmann constant. The SDF 
website provides the photospheric radii directl>3 and connects T e ff to magnitudes using 
Kenyon and Hartmann's (1995) Table A5. To present the Y2 results with (V — K) as the 
independent variable, we interpolated their T e ff to the Kenyon and Hartmann scale. The 
diameter measurements indicate, for both stars, ages in the range 10-15 MY and masses 1.6 
to 1.8 M Q . The left-hand portion of Table 3 lists more precise values obtained by placing 
the observed diameters on a finer mesh in the $ v s (V — K) color diagram (henceforth the 
$CD); the values lie in the columns designated &sdf and $y 2 . The age difference between 
HIP 560 and 21547 is not statistically significant and the SDF and Y2 tracks yield values 
that are in good agreement. The average age of the two stars is 13 ± 2 MY. The masses 
determined using the SDF and Y2 tracks are also consistent. The average mass of HIP 560 
is 1.65 ± 0.02 M Q and that of HIP 21547, 1.75 ± 0.05 M . 

If we did not have angular diameter measurement of HIP 560 and 21547, the only way 
to estimate their age and mass would be by their positions on an HRD relative theoretical 



3 www.astro.yale.edu/demarque/yystar.html 
4 www-astro. ulb.ac.be/~siess/ 
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isochrones. Figure 4 shows such an analysis, here on an HRD plotted as M K vs (V — K). 
It is seen that both the SDF and Y2 tracks indicate ages at least 5 MY older, and masses 
~ 0.2 M smaller than those indicated by the <&CD (the columns designated HRDs/^f and 
HRD Y2 in Table 4). 

The observational inputs, angular diameters, distances, and photometry, and hence 
the stellar radii, luminosities, and effective temperatures give essentially the same ages 
and masses whether the SDF or Y2 models are used. The results are however consistently 
discrepant when interpreted on the QCD or the HRD. This suggests that the source of the 
discrepancy lies in the application of theoretical calculations of non-rotating stars to stars 
that rotate with high angular velocities. 

3.2. The Effects of Stellar Rotation 

The rotation of stars less massive than ~ 2 M Q slows over their lifetimes because 
stellar winds powered by the convective outer layers carry away their angular momentum. 
The convective zone in main-sequence F spectral type stars disappears toward the earlier 
F spectral type sub-classes and energy transport becomes entirely radiative. This explains 
the rapid decrease of vsin i from FOV to F9V (e.g. Abt and Hunter 1962; Kraft 1967). The 
observation and theoretical analysis of stellar rotation have a rich history (e.g. Tassoul 2000) 
and its effects are revealed beautifully now that stars can be imaged interferometrically 
(e.g. Peterson et al. 2006 a,b; Monnier et al. 2007; Zhao et al. 2009; Che et al. 2011). 

The effective surface gravity, g e ff in a rotating star decreases from the pole to the 
equator. This produces oblateness and a brightness variation with latitude known as gravity 
darkening. The oblateness of a rotating star in radiative equilibrium is given by 
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o = 



Re 



= 0.77 




where R e and R p are the equatorial and polar radii, u is the angular velocity, and p m is the 
mean density (von Zeipel 1924; Chandrasekhar 1933). Inserting numerical values, 



where R p is the polar radius of the star, and M* and V eq are its mass and equatorial 
velocity. Sackmann (1970) showed that the decrease of R p from its non-rotating value 
is only a few percent even when the star is rotating nearly at break-up. At break-up, 
R eq = 3/2R p ,o = 0.33 and V? qbk = 2GM*/3R P (eg Collins, 1963). It is safe to apply these 
results to HIP 560 and 21547 because Demarque and Roeder (1967) showed that the outer 
convective layers disappear at F early spectral type. 

To make numerical estimates of the expected oblateness for HIP 560 and 21547, we 
use values from SDF for the radius and mass of an early F spectral type star at age 13MY, 
R* = 1.7R.0, M* = 1.7M . We use these values because, in principle, our measured values 
could be affected by oblateness. We also use the vsin i values, 171 and 95 km/s, measured 
for HIP 560 and 21547, respectively (Table 1). Since the inclinations are not known we 
cannot calculate the equatorial velocities. We can, however, estimate them by calculating 
the expectation values (vsin i). The upper bound on the inclination is i — ir/2. If the lower 
bound on i were 0, {vsin i) = V e9 7r/4. Here, however, a lower bound, icrit, is set by the 
value at which vsin i crit would imply rotation at the breakup velocity. For the adopted 
radius and mass, V eqt bk = 357 km/s. Hence, sin i crit = (observed v sini)/357 and i crit = 8.6° 
and 15.4° for HIP 560 and 21547 respectively. Using these values as a lower bound on 
the inclination, the expectation values of the equatorial velocities are 202 and 118 km/s, 
and of the oblateness, 0.104 and 0.036, for HIP 560 and 21547 respectively. These values 



o = 0.0261(^)(^)( 

V R A M/ v 100£;m/s 
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could be measured only if the stars' rotation axes were positioned for best resolution of 
the oblateness with the east-west baseline of the CHARA array. The measured oblateness 
values are therefore likely to be smaller. In the case of HIP 21547 our observations probably 
would not have measured o < 0.036 because the precision of the diameter measurement 
is ~ 2% (Table 3). For HIP 560, in the unlikely case that a ~ 10% oblateness accounts for 
the measured value of its angular diameter, a 10% smaller value of $ would still place it 
between the 10 and 20 MY isochrones. We conclude that stellar oblateness probably does 
not enter into the interpretation of our result jj. 

In a rotating star in radiative equilibrium, the radiative flux at a point on the 
photosphere is proportional to g> e // (von Zeipel 1924). Since g e ff decreases from the pole to 
the equator, the corresponding decrease in radiative flux can be characterized as a decrease 
in the local effective temperature and is called gravity darkening. In the A spectral type 
stars that have been imaged (Peterson et al. 2006a,b; Monnier et al. 2007; Zhao et al. 
2009), the temperature difference is large, ~ 2000K. Thus, relative to a non-rotating star, 
a rotating star seen pole-on (i = 0) will appear brighter and its photometric color will be 
hotter than one seen equator-on (i = vr/2). Observed values for a given star depend, of 
course, on its inclination. Maeder and Peytremann (1970) have calculated the photometric 
effects for stars of mass 1.4, 2, and 5 M Q at various values of V eq and i. Typically the values 
of My vary by a few tenths mag and of (B — V) by a few hundreths. Details depend on the 
stellar mass and inclination; Maeder and Peytremann's calculations show that for i > 54° 
the rotating star appears less bright and redder than a non-rotating star. 

The effects of gravity darkening suggest that the discrepancies summarized in Table 
4 are attributable to rotation of HIP 560 and 21547. If their rotational velocities and 



5 Also, the on-sky angle of the projected baseline was essentially the same throughout our 
observations (Table 2). Thus baseline rotation during our observations is not an issue. 
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inclinations are such that their absolute magnitudes at K are depressed by a few tenths 
relative to a non-rotating star and their (V — K) are slightly redder, the ages and masses 
deduced from the <&CD and HRD could be in agreement. A detailed test of this hypothesis 
will be possible when the inclinations of these stars are measured. This can be accomplished 
by either measuring their light fluctuations and thus rotational periods or by interferometric 
imaging of their photospheric emission of these stars is mapped. The first approach is 
possible now as Garcia-Alvarez et al. (2011) demonstrated by measuring the periods and 
inclinations of two other stars in the BPMG. 

3.3. Implications for the BPMG 

If we accept the results of the <&CD analysis and attribute the discrepancy of results 
from the HRD to stellar rotation, we conclude that HIP 560 and 21547 are roughly coeval 
at an average age of about 13 MY. Their galactic (X,Y,Z) coordinates (Table 1) place them 
about 49 pc apart with each moving at speed ~ 22 km/s given by their (U,V,W) velocity 
components. It is too soon to attribute the results for this pair to all the stars in the BPMG 
but our results suggest that the members were born together about 13 MY ago, and thus 
are at the younger end of the age span described in §2. 

4. Summary and Future Directions 

1) The measured angular diameters of HIP 560 and HIP 21547 in the H and K bands are 
0.492 ± 0.032 and 0.518 ± 0.009 mas, respectively. Scaled to a common distance lOpc, their 
angular diameters are $ = 1.94 ± 0.13 and 1.52 ± 0.03 mas. 

2) Analyzing these results in $ vs (V — K) diagram with SDF and Y2 isochrones calculated 
yields ages of the two stars that are not different at a statistically significant level. Their 
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average age is 13 ± 2 MY. The masses determined according to the two sets of isochrones 
also do not differ at a statistically different level. The average values are 1.65 ± 0.02 M for 
HIP 560 and 1.75 ± 0.05 M for HIP 21547. 

3) Analyzing the stellar parameters with the SDF and Y2 isochrones in a conventional HR 
diagram yields ages at least 5 MY older and masses ~ 0.2 M Q smaller. 

4) HIP 560 and 21547 are rapid rotators (Table 1). The discrepancy in ages and masses can 
be accounted for by gravitational darkening of rotating stars in radiative equilibrium. A 
detailed test of this hypothesis will be possible when their their inclinations are determined 
either by measuring their rotational periods or by interometric images of their photospheres. 

5) Taken together, our results suggest that HIP 560 and 21547 formed coevally about 13 
MY ago. 

The most pressing task is to determine whether the 13 MY age applies to the BPMG 
as a whole. Although most BPMG are in the southern hemisphere, at least 3 other F and 
G spectral type stars are bright and near enough for measurement with the CHARA array. 
Stars of later spectral type are too faint for observation at the present time unless, by good 
luck, they are close to the sun. Many more of the BPMG members will be observable when 
the planned sensitivity improvements at the CHARA array are realized. 

We thank the referee for a helpful and unusually warm report. We are grateful to 
Steve Ridgway for advice and help with the observations. We thank Deane Peterson for 
reminding us that F-stars can rotate rapidly, for advice about stellar rotation, and for 
pointing out the Garcia-Alvarez paper. We thank the CHARA staff at Mt. Wilson for 
thorough support. The CHARA Array is operated by Georgia State University's Center for 
High Angular Resolution Astronomy on Mount Wilson, California. Access to the CHARA 
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Table 1. Properties of HIP 560 and HIP 21547 



HIP 560 



HIP 21547 



Alternate Name 
Spectral Type 
V(mag) 
K(mag) 
vsini (km/s) 
Distance (pc) 
X,Y,Z (pc) 
U,V,W (km/s) 



HD 203 
F3V 
6.19 
5.24 
170.7 
39.4 ±0.4 
4.5, 5.8, -38.4 
-10.4, -14.5, -13.3 



HD 29391 
FOV 
5.22 
4.54 
95.0 
29.4 ±0.3 
-24.3, -8.2, -15.2 
-14.0, -16.2, -10.1 



Calibrators: 



4>LD (mas) 



(mas) 



HD 268 
0.27 
HD 223884 
0.26 



HD 26794 
0.25 



Data for /? Pic targets from Torres et al. (2006) 
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Table 2. Calibrated Visibilities 



Object 


\ 

A 


MJD 


a (m) 


PA (°\ 


Veal 


a Vcal 


HIP 


560 


K 


55450.286 


308.62 


69.69 


0.828 


0.018 


HIP 


560 


K 


55450.300 


311.91 


71.87 


0.783 


0.016 


HIP 


560 


K 


55450.314 


313.44 


73.94 


0.751 


0.015 


HIP 


560 


K 


55450.332 


312.66 


76.38 


0.866 


0.017 


HIP 


560 


K 


55450.347 


309.38 


78.47 


0.863 


0.016 


HIP 


560 


K 


55452.314 


313.51 


74.74 


0.839 


0.015 


HIP 


560 


K 


55452.324 


312.91 


76.10 


0.814 


0.015 


HIP 


560 


K 


55452.334 


311.39 


77.39 


0.856 


0.015 


HIP 


560 


H 


55451.328 


312.81 


76.22 


0.853 


0.028 


HIP 


560 


H 


55451.342 


310.18 


78.08 


0.897 


0.025 


HIP 


560 


H 


55451.355 


305.80 


79.82 


0.893 


0.025 


HIP 


560 


H 


55451.370 


298.79 


81.72 


0.866 


0.026 


HIP 


21547 


K 


55450.515 


313.26 


75.68 


0.849 


0.013 


HIP 


21547 


K 


55450.525 


313.44 


75.86 


0.819 


0.012 


HIP 


21547 


K 


55450.535 


312.45 


75.96 


0.860 


0.014 


HIP 


21547 


H 


55451.501 


311.41 


75.40 


0.769 


0.021 


HIP 


21547 


H 


55451.515 


313.41 


75.73 


0.776 


0.022 


HIP 


21547 


H 


55451.533 


312.20 


75.97 


0.766 


0.021 


HIP 


21547 


H 


55451.544 


310.15 


76.02 


0.775 


0.021 


HIP 


21547 


11 


55452.496 


311.11 


75.37 


0.773 


0.022 


HIP 


21547 


H 


55452.513 


313.43 


75.73 


0.719 


0.020 


HIP 


21547 


H 


55452.521 


313.31 


75.88 


0.791 


0.022 


HIP 


21547 


H 


55452.533 


311.69 


75.99 


0.729 


0.020 



Table 3. Measured Sizes 



HIP 560 HIP 21547 



<t>{mas) 0.492 ± 0.032 0.518 ± 0.009 
<3>(mas) 1.94 ±0.13 1.52 ± 0.03 
R/R & 2.08 ±0.14 1.63 ±0.03 
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Table 4. Measured Ages and Masses 



Star 


®SDF 


<S> Y2 


HRD 5DF 


HRDy 2 






Age (MY) 




HIP 560 


12.0 ±3.0 


10 ±2 


25±f 


18 ±2 


HIP 21547 


15.0 ±2.0 


15 ±2 


2012° 


18 ±2 






Mass 


(M ) 




HIP 560 


1.68 ± 0.03 


1.63 ±0.02 


1.49 ±0.05 


1.45 ±0.10 


HIP 21547 


1.75 ± 0.05 


1.75 ±0.05 


1.56 ±0.05 


1.55 ±0.10 
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Fig. 1. — Calibrated visibilities of HIP 560 at H and K vs spatial frequency, the projected 
baseline divided by the wavelength. The dashed curves are for uniform disk models (UD) 
with the diameters indicated and the solid curves are limb-darkened models. 
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Fig. 2.— Same as Fig. 1 but for HIP 21547. 
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Fig. 3. — HIP 560 and 21547 angular diameters scaled to lOpc distance vs (V-K) compared 
to isochrones at 10, 20, and 50 MY calculated by SDF (top) and Y2 (bottom). The crosses 
indicate masses 0.4 to 2.0 M at 0.1 M intervals for the SDF isochrones and 0.5 to 2.0 
M Q for the Y2 isochrones. 



-22 - 




8 I , , , I , , , I , , , I 

2 4 6 

V-K 

Fig. 4. — HIP 560 and 21547 on a conventional Hetrzsprung-Russell diagram here presented 
as Mk vs (V-K). The isochrones are calculated by SDS and Y2 as in Fig. 3. 



